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We study charm production in ultra-relativistic heavy-ion collisions by using the Parton-Hadron- 
String Dynamics (PHSD) transport approach. The initial charm quarks are produced by the Pythia 
event generator tuned to fit the transverse momentum spectrum and rapidity distribution of charm 
quarks from Fixed-Order Next-to-Leading Logarithm (FONLL) calculations. The produced charm 
quarks scatter in the quark-gluon plasma (QGP) with the off-shell partons whose masses and widths 
are given by the Dynamical Quasi-Particle Model (DQPM), which reproduces the lattice QCD 
equation-of-state in thermal equilibrium. The relevant cross sections are calculated in a consistent 
way by employing the effective propagators and couplings from the DQPM. Close to the critical 
energy density of the phase transition, the charm quarks are hadronized into D mesons through 
coalescence and fragmentation. The hadronized D mesons then interact with the various hadrons 
in the hadronic phase with cross sections calculated in an effective lagrangian approach with heavy- 
quark spin symmetry. Finally, the nuclear modification factor Raa and the elliptic flow V 2 of D° 
mesons from PHSD are compared with the experimental data from the STAR Collaboration for 
Au+Au collisions at ^/snn =200 GeV. We find that in the PHSD the energy loss of D mesons at 
high pr can be dominantly attributed to partonic scattering while the actual shape of Raa versus 
Pt reflects the heavy-quark hadronization scenario, i.e. coalescence versus fragmentation. Also the 
hadronic rescattering is important for the Raa at low pr and enhances the D-meson elliptic flow 
V2- 
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I. INTRODUCTION 

According to the fundamental theory of stronginterac- 
tions, the Quantum Chromo Dynamics (QCD) [jj, mat¬ 
ter changes its phase at high temperature and density and 
bound (colorless) hadrons dissolve to interacting (col¬ 
ored) quarks and gluons - Quark-Gluon-Plasma (QGP). 
Such extreme conditions have existed in the early ex¬ 
pansion of the universe and now can be realized in the 
laboratory by collisions of heavy-ions at ultra-relativistic 
energies. The study of the phase boundary and the prop¬ 
erties of the QGP are the main goal of several present and 
future heavy-ion experiments at SPS (Super Proton Syn¬ 
chrotron), RHIC (Relativistic Heavy-Ion Collider), LHC 
(Large Hadron Collider) and the future FAIR (Facility 
for Antiproton and Ion Research) and NICA (Nuclotron- 
based Ion Collider fAcility) [2J. Since the QGP is cre¬ 
ated only for a short time (of a couple of fm/c) it is 
quite challenging to study its properties and to find the 
most sensible probes. In order to study the full com¬ 
plexity of the underling problem, one needs to obtain 
comprehensive information by the measurement of the 
’bulk’ light hadrons, electromagnetic probes (dileptons 
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and photons), heavy mesons and jets. The advantage of 
the ’hard probes’ such as the mesons containing heavy 
quarks (charm and beauty) is, firstly, that due to the 
heavy masses they are dominantly produced in the very 
early stages of the reactions with large energy-momentum 
transfer, contrary to the light hadrons and electromag¬ 
netic probes. Secondly, they are not in an equilibrium 
with the surrounding matter due to smaller interaction 
cross sections relative to the light quarks and, thus, 
may provide an information on their creation mecha¬ 
nisms. Moreover, due to the hard scale, perturbative 
QCD (pQCD) should be applicable for the calculation 
of heavy quark production. As shown in Ref. Q, the 
FONLL calculations are in good agreement with the ex¬ 
perimental observables on charm meson spectra in p+p 
collisions. This provides a solid reference frame for study¬ 
ing the heavy-meson production in heavy-ion collisions. 

The collective properties of open charm mesons have 
been addressed experimentally by measuring the nuclear 
modification factor Raa, which is the ratio of the trans¬ 
verse momentum distribution in A + A collisions relative 
to p+p collisions scaled by the number of binary colli¬ 
sions, as well as the collective elliptic flow V 2 ■ These 
two observables are comprehensive since the high pp part 
of Raa is very sensitive to the energy loss of charm 
quarks (or mesons) during their propagation through 
the partonic (or hadronic) medium due to the interac- 
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tion processes. The low px part is more sensitive to the 
hadronization mechanism and, thus, provides constraints 
on the relative scale for the transition from partonic to 
hadronic degrees-of-freedom. Moreover, the elliptic flow 
V 2 of charm quarks characterizes the collectivity of the 
system developed from the very early stage to the freeze- 
out time thereby including interactions from the partonic 
and hadronic phase. Thus, the final momentum distribu¬ 
tion of the heavy mesons is sensitive to the details of the 
expansion of the plasma itself as well as to the strength of 
the interaction of the heavy quark with the partons and 
the formed D-mesons with the hadronic environment. It 
depends on the density of partons (and their properties) 
present in the QGP and hence on whether the plasma is 
in equilibrium during the expansion. 

The first charm measurements at RHIC energies by the 
PHENIX0 and STAR @ collaborations were related to 
the single non-photonic electrons emitted from the decay 
of charm mesons. However, recently the STAR Collabo¬ 
ration measured directly the nuclear modification factor 
and the elliptic flow of D° mesons in Au+Au collisions at 
^/snn =200 GeV 0, ;7] which allows for a straight forward 
comparison with the theoretical model calculations. It 
has been observed that the Raa and V 2 of charm mesons 
show a similar behavior as in case of light hadrons con¬ 
trary to expectations from pQCD. Similar observations 
have been made at LHC energies, too 

It still remains a challenge for the theory to reproduce 
the experimental data and to explain simultaneously the 
large energy loss of charm quarks (Raa) and the strong 
collectivity (i^) - cf. e.g. Refs. |9l-|25|. The interactions 
of charm quarks with the partonic medium are commonly 
based on pQCD with massless light quarks and a fixed 
or running coupling. The time evolution of the charm- 
quark distribution in the expanding fireball is approxi¬ 
mated by the Fokker-Plank equation where the response 
of the partonic (or hadronic) medium is expressed in 
terms of temperature- and momentum-dependent drag 
and diffusion coefficients. The modeling of the tempera¬ 
ture profile is often done in a fireball model or by hydro- 
dynamic calculations (ideal or viscous) which start with 
some initial conditions and follow the dynamical evolu¬ 
tion according to the chosen Equation-of-State (EoS) un¬ 
der the assumption of local equilibrium. In partonic cas¬ 
cade models one solves the Boltzmann equation for mass¬ 
less quarks/gluons with the pQCD cross sections for some 
fixed coupling a s . The hadronization of charm quarks is 
done assuming coalescence at low px and fragmentation 
at high pt- In spite that many models may describe the 
Raa, it is still difficult to obtain simultaneously a con¬ 
sistent description of the elliptic flo w vo. using the same 
assumptions and model parameters [251 1 - e -S- the choice 
of the running coupling, the K-factor to scale the pQCD 
cross sections, the time evolution profile, etc. Moreover, 
the conclusions on the amount of suppression due to col- 
lisional energy loss by means of the elastic interactions of 
charm quarks with the QGP partons versus the radiative 
energy loss due to the emission of soft gluons (i.e. gluon 


bremsstrahlung) are still far from being robust. Also the 
influence of hadronization and especially hadronic rescat¬ 
tering is not yet settled, too. Moreover, the results turned 
out to be also sensitive to the time evolution models in¬ 
volved. 

Our goal here is to study the charm dynamics based 
on a consistent microscopic transport approach for the 
charm production, hadronization and rescattering with 
the partonic and hadronic medium. In this study we 
will confront our calculations within the Parton-Hadron- 
String Dynamics (PHSD) approach to the experimental 
data on charm at RHIC energies and discuss the per¬ 
spectives/problems of using the charm quarks for the 
tomography of the QGP. To achieve this goal we em¬ 
bed the heavy-quark physics in the existing PHSD trans¬ 
port approach 1261 which incorporates explicit partonic 
degrees-of-freedom in terms of strongly interacting quasi¬ 
particles (quarks and gluons) in line with an equation- 
of-state from lattice QCD (1QCD) as well as dynamical 
hadronization and hadronic elastic and inelastic collisions 
in the final reaction phase. Since PHSD has been success¬ 
fully applied to describe the final distribution of mesons 
(with light quark content) from lower SPS up to LHC 
energies [20 - 1291 ] , it provides a solid framework for the de¬ 
scription of the creation, expansion and hadronization of 
the QGP as well as the hadronic expansion with which 
the heavy quarks interact either as quarks or as bound 
states such as D-mesons. 

Thus, the principle differences of our approach to the 
previous dynamical models (including earlier HSD [30) 
studies on the charm dynamics [3lJ) are: 

(i) the degrees-of-freedom for the QGP are massive and 
strongly interacting quasiparticles contrary to massless 
and weakly interacting pQCD partons; 

(ii) a non-equilibrium off-shell microscopic transport 
approach is employed for the QGP dynamics, hadroniza¬ 
tion and the hadronic phase instead of simplified 
descriptions of the parton dynamics in terms of Fokker- 
Planck equations + QGP hydrodynamics (assuming 
local equilibration) or Boltzmann-type partonic cascades 
with massless light quarks; 

(iii) rescattering of D-mesons in the hadronic phase in 
line with an up-to-date effective Lagrangian approach 
from Ref. [32| . 

This paper is organized as follows: In Sec. [IT] we recall 
the basic ideas of the PHSD approach while in Sec. |m| 
we describe the production of initial charm quark pairs 
in hard binary nucleon-nucleon collisions and their im¬ 
plementation in PHSD. In Sec. m we present the in¬ 
teractions of charm quarks with off-shell partons in the 
QGP calculated earlier in Refs. [3314351 ] . the hadroniza¬ 
tion of charm quarks in Sec. El and the D meson inter¬ 
actions with hadrons in Sec. IVII that are based on the 
cross sections from Ref. (32|. Finally, the nuclear modi- 
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fication factor Raa and the elliptic flow V 2 of D mesons 
from PHSD are presented in Sec. IVIII and compared to 
the available data. A summary completes this work in 

Sec. Enu 


II. THE PHSD TRANSPORT APPROACH 

The Parton-Hadron-String Dynamics (PHSD) trans¬ 
port approach [2(1 [H} is a microscopic covariant dy¬ 
namical model for strongly interacting systems formu¬ 
lated on the basis of Kadanoff-Baym equations [HI [33] 
for Green’s functions in phase-space representation (in 
first order gradient expansion beyond the quasiparticle 
approximation). The approach consistently describes 
the full evolution of a relativistic heavy-ion collision 
from the initial hard scatterings and string formation 
through the dynamical deconfinement phase transition to 
the strongly-interacting quark-gluon plasma (sQGP) as 
well as hadronization and the subsequent interactions in 
the expanding hadronic phase as in the Hadron-String- 
Dynamics (HSD) transport approach [3(3]. The trans¬ 
port theoretical description of quarks and gluons in the 
PHSD is based on the Dynamical Quasi-Particle Model 
(DQPM) for partons that is constructed to reproduce 
1QCD results for a quark-gluon plasma in thermodynamic 
equilibrium [HI on the basis of effective propagators for 
quarks and gluons. The DQPM is thermodynamically 
consistent and the effective parton propagators incorpo¬ 
rate finite masses (scalar mean-fields) for gluons/quarks 
as well as a finite width that describes the medium de¬ 
pendent reaction rate. For fixed thermodynamic parame¬ 
ters (T, p q ) the partonic width’s ri(T, p q ) fix the effective 
two-body interactions that are presently implemented in 
the PHSD [13]. The PHSD differs from conventional 
Boltzmann approaches in a couple of essential aspects: 
i) it incorporates dynamical quasi-particles due to the fi¬ 
nite width of the spectral functions (imaginary part of the 
propagators); ii) it involves scalar mean-fields that sub¬ 
stantially drive the collective flow in the partonic phase; 
iii) it is based on a realistic equation of state from lattice 
QCD and thus describes the speed of sound c s (T) reli¬ 
ably; iv) the hadronization is described by the fusion of 
off-shell partons to off-shell hadronic states (resonances 
or strings) and does not violate the second law of thermo¬ 
dynamics; v) all conservation laws (energy-momentum, 
flavor currents etc.) are fulfilled in the hadronization 
contrary to coalescence models; vi) the effective partonic 
cross sections no longer are given by pQCD and are ’de¬ 
fined’ by the DQPM in a consistent fashion and probed 
by transport coefficients (correlators) in thermodynamic 
equilibrium by performing PHSD calculations in a fi¬ 
nite box with periodic boundary conditions (shear- and 
bulk viscosity, electric conductivity, magnetic suscepti¬ 
bility etc. 0,53). 

In the beginning of relativistic heavy-ion collisions 
color-neutral strings (described by the FRITIOF LUND 
model [13]) are produced in hard scatterings of nucle¬ 


ons from the impinging nuclei. These strings are dis¬ 
solved into ’pre-hadrons’ with a formation time of ~ 0.8 
fm/c in their rest frame, except for the ’leading hadrons’, 
i.e. the fastest residues of the string ends, which can re¬ 
interact (practically instantly) with hadrons with a re¬ 
duced cross sections in line with quark counting rules. 
If, however, the local energy density is larger than the 
critical value for the phase transition, which is taken to 
be ~ 0.5 GeV/fm 3 , the pre-hadrons melt into (colored) 
effective quarks and antiquarks in their self generated re¬ 
pulsive mean-field as defined by the DQPM [HI- In the 
DQPM the quarks, antiquarks and gluons are dressed 
quasiparticles and have temperature-dependent effective 
masses and widths which have been fitted to lattice ther¬ 
mal quantities such as energy density, pressure and en¬ 
tropy density. The nonzero width of the quasiparticles 
implies the off-shellness of partons, which is taken into 
account in the scattering and propagation of partons in 
the QGP on the same footing (i.e. propagators and cou¬ 
plings). We point out that the DQPM does not provide 
effective propagators for the c— and c—quarks since the 
latter degrees-of-freedom are subdominant in the entropy 
density for temperatures of 1 - 3 T c due to their large 
mass and thus cannot be determined properly by the en¬ 
tropy from the lattice. Nevertheless, the interactions of 
the charm quarks with the light (it, d, s) quarks can be 
computed on the basis of the running coupling g 2 (T/T c ) 
in the DQPM and the effective pro pag otors for the light 
quarks. As demonstrated in Ref. (33 the resulting dif¬ 
ferential cross sections of c, c quarks only very slightly 
depend on the spectral width of the charm quarks such 
that even the on-shell limit for these degrees-of-freedom 
with a mass of ~ 1.5 GeV provides a very reasonable 
approximation [HI HU. 

The transition from the partonic to hadronic degrees- 
of-freedom (for light quarks/antiquarks) is described 
by covariant transition rates for the fusion of quark- 
antiquark pairs to mesonic resonances or three quarks 
(antiquarks) to baryonic states, i.e. by the dynami¬ 
cal hadronization. We already mention here that this 
hadronization process is restricted to ’bulk’ transverse 
momenta pr up to ~ 2 GeV and has to be replaced by 
fragmentation for high pt . Note that due to the off-shell 
nature of both partons and hadrons, the hadronization 
process described above obeys all conservation laws (i.e. 
four-momentum conservation and flavor current conser¬ 
vation) in each event, the detailed balance relations and 
the increase in the total entropy S. In the hadronic 
phase PHSD is equivalent to the hadron-strings dynam¬ 
ics (HSD) model [30] that has been employed in the 
past from Schwerlonen-Synchrotron (SIS) to SPS ener¬ 
gies. On the other hand the PHSD approach has been 
applied to p+p, p+A and relativistic heavy-ion collisions 
from lower SPS to LHC energies and been successful in 
describing a large number of experimental data including 
single-particle spectra, collective flow as well as electro¬ 
magnetic probes [26l-[29f. 

In this study we use the Pythia event generator [43} 
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for nucleon-nucleon binary collisions to produce charm- 
quark pairs in relativistic heavy-ion collisions and evolve 
the dynamics of the charm quarks within PHSD in or¬ 
der to understand the dominant mechanisms in compar¬ 
ison to the recent experimental data on D° mesons from 
STAR. Whereas the dynamics of the charm quarks is 
essentially Boltzmann-like, i.e. without mean-fields and 
dynamical width, their interactions with the dynamical 
partons (light quarks and gluons) is based on the effective 
coupling from the DQPM and the DQPM propagators 

sip. 


III. INITIAL CHARM QUARK PRODUCTION 

Before studying the charm production in relativistic 
heavy-ion collisions, we discuss the charm production in 
p+p collisions at the top RHIC energy of ^/snn = 200 
GeV. The charm production in p+p collisions also plays 
the role as a reference for the nuclear modification factor 
Raa in heavy-ion collisions. 

We use the Pythia event generator to produce charm 
and anticharm quarks in p+p collisions with the pa¬ 
rameters PARP(91)=1.0 GeV/c and PARP(67)=1.0 as 
in Ref. [S [p • The former parameter is the Gaussian 
width of the primordial transverse momentum of a par- 
ton which initiates a shower in hadrons and the latter 
the parton shower level parameter. We note that by 
an additional suppression of the transverse momenta of 
charm and anticharm quarks by 10 % and a suppression 
of rapidities by 16 % the transverse momentum spectrum 
as well as rapidity distribution of charm and anticharm 
quarks from the Pythia event generator are very similar 
to those from the FONLL calculations in p+p collisions 
at -y/sNN =200 GeV as shown in Fig. |T] Here the red 
dot-dashed lines are from the Pythia event generator (af¬ 
ter tuning) and the blue dotted lines from the FONLL 
calculations, respectively. 

The produced charm and anticharm quarks hadronize 
by emitting soft gluons. The probabilities for a charm 
quark to hadronize into D°, D + , D *°, D* + , Df, and 
A c are, respectively, taken to be 0.2, 0.174, 0.213, 0.224, 
0.08, and 0.094 (45|. The momentum of the hadronized 
D meson or A c is given by the fragmentation func¬ 
tion |46|, 
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FIG. 1: (Color online) Transverse momentum (upper) and ra¬ 
pidity (lower) distributions of charm quarks in p+p collisions 
at -y/sNN =200 GeV from FONLL (dotted lines) and the tuned 
Pythia event generator (dot-dashed lines); in the upper part 
the transverse momentum spectrum of D° mesons which are 
fragmented from charm quarks with the contribution from D * 
decay included (solid line) and that of D* + after scaling are 
compared with the experimental data from the STAR Collab¬ 
oration fill ]. 
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where z is the momentum fraction of the hadron H frag¬ 
mented from the heavy quark Q while cq is a fitting pa¬ 
rameter which is taken to be £q = 0.01 in our study. We 
note that the parameter cq used here is smaller than the 
usual value because the transverse momentum of charm 
quarks is reduced in our study. The solid line in Fig. [T] 
(upper part) shows the transverse momentum spectrum 
of D° mesons after charm quark fragmentation including 
the contribution from the decay of D*° and D* + . We 


can see that our results reproduce the experimental data 
from the STAR Collaboration [44} reasonably well. 

In heavy-ion collisions, the nucleons in the target and 
projectile nuclei are distributed in coordinate space ac¬ 
cording to the nuclear density distribution given by 


1 + exp[(r - r 0 )/a\ ’ 

where r 0 = 1.124A 1 / 3 and a = 0.02444A 1 / 3 + 0.2864 with 
A being the mass number of the target or projectile nu¬ 
cleus. Each nucleon has Fermi motion depending on the 
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local nucleon density which is chosen randomly by Monte 
Carlo. Although the Fermi momentum is small in the rest 
frame of each nucleus, the component along the beam 
direction becomes large in the laboratory frame due to 
the Lorentz transformation ( 7 cm « 100 at the top RHIC 
energy). Accordingly, the Fermi motions smears the en¬ 
ergy of nucleon-nucleon binary collisions in relativistic 
heavy-ion collisions. The lower part of Fig. [2] shows the 
distribution of binary nucleon-nucleon collision energies 
in 0-10 % central Au+Au collisions at ^/snn =200 GeV 
from PHSD. We mention that this Fermi smearing is usu¬ 
ally neglected in theoretical models for charm production 
and propagation which is not crucial at RHIC energies 
but becomes important close to threshold energies for 
charm production. 

Since charm-quark production requires a high energy- 
momentum transfer, the number of produced charm 
quark pairs in relativistic heavy-ion collisions is propor¬ 
tional to the number of binary nucleon-nucleon collisions 
Nbi n . Since the probability to produce a charm quark 
pair depends on invariant energy and is less than that for 
primary hard collision in the PHSD, the binary nucleon- 
nucleon collisions producing charm quark pairs are cho¬ 
sen by Monte Carlo from the ratio of the cross section for 
charm production in nucleon-nucleon collision, a^g(y/s), 
to the inelastic nucleon-nucleon cross section. The total 
cross section in PHSD for charm production in nucleon- 
nucleon collisions is parameterized by 

= A[1 - >/(V^)]“(V* of, (3) 

where A, a and /3 are fit parameters. We use two different 
parameter sets, one for y/s <200 GeV and the other for 
y/s >200 GeV to fit the experimental data covering a 
wide range of collision energies as shown in the upper part 
of Fig. [2] both parameterizations are smoothly connected 
at y/s =200 GeV. 

The energy-momenta of the produced charm and an¬ 
ticharm quarks in each collision event are given by the 
Pythia event generator as in p+p collisions but for the 
actual (smeared) collision energy. In the Pythia event 
generator, the charm and anticharm quarks are produced 
in the center-of-mass frame of two colliding nucleons with 
the incident nucleons being aligned along the z-direction. 
Therefore in PHSD we rotate the momenta of the gener¬ 
ated charm and anticharm quark to the original orienta¬ 
tion of the two incident nucleons in their center-of-mass 
frame and then boost to the calculational frame. We as¬ 
sume that the two nucleons, which produce the charm 
quark pair, keep their transverse momentum and loose 
only longitudinal momentum such that the total energy 
(including the produced charm quark pair) is conserved 
and the two nucleons are still on-shell after the hard colli¬ 
sion. Although this prescription seems to violate spatial 
momentum conservation, we have to mention that this 
violation is very small since the charm quark pair is pro¬ 
duced together with several (plenty) light hadrons in the 
same event which balance the spatial momentum. 




FIG. 2: (Color online) The total cross section for charm pro¬ 
duction in p+p collisions (as parameterized in the PHSD) is 
compared with the experimental data at various collision ener¬ 
gies I44l. l47l (upper panel); the distribution of binary nucleon- 
nucleon collision energies in 0-10 % central Au+Au collisions 
at y /snn =200 GeV from PHSD including Fermi smearing 
(lower panel). 


IV. CHARM-QUARK SCATTERING IN THE 
QGP 

A. q + c and g + c elastic scattering cross sections 

Quarks, antiquarks, and gluons are dressed in the QGP 
and have temperature-dependent effective masses and 
widths. In the DQPM, the mass and width of the light 
partons are given by thermal quantum-field theory as¬ 
suming leading order diagrams but the strong coupling 
g 2 (T) is fitted to lattice data on energy and entropy den¬ 
sities, etc. [38]. Note that a nonzero width of a parton 
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reflects the off-shell nature as well as the strong interac¬ 
tion of the quasi-particle and finite life-time. 

The charm and anticharm quarks produced in early 
hard collisions interact with the dressed off-shell partons 
in the QGP. The cross sections for the charm quark scat¬ 
tering with massive off-shell partons are calculated con¬ 
sidering the mass spectra of final state particles 0, 0 . 
In the current study the charm quark mass is taken to be 
1.5 GeV and its mass spectrum is neglected for simplic¬ 
ity. Since the charm quark is heavy, the off-shell effect is 
small contrary to the light quarks and gluons (see Refs. 
0 0 for a quantitative study of this effect). More¬ 
over, the off-shell effects from parton spectral functions 
are small, except for reducing the kinematic threshold. 



In the current study we consider only elastic scatter¬ 
ing of charm quarks by light quarks and gluons. We do 
not consider yet the radiative processes which generate 
radiative energy loss because we expect that, due to the 
large gluon mass in the DQPM, the radiative processes 
are sub-dominant as compared to the collisional ones, es¬ 
pecially for low charm-quark momenta (pr ). We expect 
the radiative energy loss to contribute at very high pp as 
accessible experimentally at the LHC [0. 

The elastic scattering of charm quarks in the QGP is 
treated in our study by including the non-perturbative 
effects of the strongly interacting quark-gluon plasma 
(sQGP) constituents, i.e. the large coupling, the mul¬ 
tiple scattering etc. The multiple strong interactions of 
quarks and gluons in the sQGP are encoded in their ef¬ 
fective propagators with broad spectral functions. As 
pointed out above the effective propagators, which can 
be interpreted as resummed propagators in a hot and 
dense QCD environment, have been extracted from lat¬ 
tice data in the scope of the DQPM [38']. Furthermore, 
in Refs. 0,0 we have evaluated, to the lowest order in 
the perturbation expansion, the transition amplitudes of 
the processes q + c —> q + c and g + c—> g + c a,s shown in 
Fig. [3] Contrary to the case of massless gluons where the 
“Transverse gauge” is used, we use the “Lorentz covari¬ 
ance” for the case of massive gluons since a finite mass in 
the gluon propagator allows to fix the O’th components 
of the gluon fields A° a (a = 1;...; 8) by the spatial degrees 
of freedom A^(k = 1; 2;3). Furthermore, the divergence 
encountered in the f-channel when calculating the total 
cross sections (<r qc ) and (<r gc ) is cured self-consistently 
in our computation since the infrared regulator is given 
by the finite DQPM gluon mass (and width). 

For the scattering of charm quark by the light quark 
and gluon with finite masses and widths, we take into ac¬ 
count the spectral functions of the light quark and gluon, 
and the temperature-dependent running coupling, g 2 (T). 
The propagators of massive vector gluons with finite life¬ 
time and of charm quark with zero life time are, respec¬ 
tively, given by 


FIG. 3: (Color online) Feynman diagrams for the scattering 
of g + c —> q + c and g + c —> g + c. Latin (Greek) subscripts 
denote colour (spin) indices, ki and pi ( kf and pf) denote 
the initial (final) 4-momentum of the light quark or the gluon 
and the heavy quark, respectively. 


G^(q) 

S(P) 


— /--- 

q 2 -m 2 g + i2"f g q 0 ’ 

if + m c 

9 9 5 

p z — m z 


(4) 


where m g , are Respectively, gluon mass and width 
at finite temperature, and the charm quark mass m c is 
taken to be 1.5 GeV. 

Fig. 0] shows the elastic cross section of an on-shcll 
charm quark with off-shell u-quarks (tr uc ) as a func¬ 
tion of the temperature and yfs, the energy available in 
the center-of-mass frame. Apart from the threshold re¬ 
gion the cross section is rather independent of y/s, and 
decreases with increasing temperature mainly due to a 
smaller coupling strength at high temperature. The large 
cross section near the critical temperature T c = 0.158 
GeV is related to the infrared enhancement of the cou¬ 
pling a s (T ) in the DQPM. The cross section for charm 
quark and gluon elastic scattering is about twice (9/4 
which is the ratio of the different color Casimir operators 
squared) that for charm quark and light quark scattering. 

Our cross sections also depend on the scattering angle 
besides collision energy yfs and temperature T. Fig. [5] 
presents the differential cross sections da/d cos 9 for on- 
shell charm quark scattering with off-shell u-quarks at 
yfs = 3 GeV (orange lines) and 4 GeV (black lines) 
for temperatures of 1.2 T c , 2 T c and 3 T c . Different 
from the pQCD-inspired models, da/dcosd in our ’non- 
perturbative’ model is not so much forward peaked (large 
enhancement for small angles or small momentum trans¬ 
fers t) [33]. Therefore, efficient momentum transfers more 
often occur in our partonic scattering as compared to the 
usual pQCD 2 —>• 2 scattering [34]. 
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FIG. 4: (Color online) Elastic cross section for u + c — > u + c 
scattering as a function of the temperature T and the invari¬ 
ant energy y/s where the M-quark is off-shell and the charm 
quark has a constant mass of 1.5 GeV [33. .341], 



FIG. 5: (Color online) Differential elastic cross section for the 
scattering of an on-shell charm quark and off-shell u-quarks at 
y/s = 3 GeV (orange lines) and 4 GeV (black lines) for tem¬ 
peratures of 1.2 T c , 2 T c and 3 T c , with T c = 0.158 GeV. We 
take 1.5 GeV for the charm quark mass and DQPM spectral 
functions (propagators) for the light off-shell partonic degrees- 
of-freedom [331 ]. 


B. Charm spatial diffusion coefficient 


Using the transition amplitudes for the elastic scatter¬ 
ing of charm quarks by the partons in medium, one can 
calculate the charm spatial diffusion coefficient D s from 


the drag coefficient, A = d< £ c> through tjd = A/p c : 


De = lim 


T 


Pc^o m c r] D 


(5) 


or from the diffusion coefficient k = i tf<(Pc „ p - 1 > : 

3 at 


d, = lim 


Pc-> o 


( 6 ) 


Both definitions ag ree with each other if the Einstein 
relation is valid [34l . l35| . 

Using Eq. (0, we show in Fig|6]the spatial diffusion co¬ 
efficient D s as a function of T for = 0. Our results are 
compared to the lattice calculations, which have recently 
been confirmed by the Bielefeld collaboration [49j], for 
temperatures above T c ss 160 MeV and with the spatial 
diffusion coefficient of a heavy meson in hadronic mat¬ 
ter [3)3] for temperatures below 180 MeV. We observe that 
the spatial diffusion coefficients in hadronic and partonic 
matters meet each other and have a pronounced mini¬ 
mum around T c . Fig.[6]shows that our results agree with 
those from the lattice QCD. The smooth transition of the 
heavy-quark transport coefficients from the hadronic to 
the partonic medium corresponds to a crossover in line 
with lattice calculations, and differs substantially from 
perturbative QCD (pQCD) calculations which show a 
large discontinuity of D s close to T c [H, Q. 

We emphasize that the transport coefficient D s ex¬ 
tracted from our microscopic calculations and its agree¬ 
ment with the 1QCD results (within errors) and the cor¬ 
responding D meson D s in hadronic medium validate 
our description for the coupling of charm with the QGP 
matter. 

Fig- 0 shows the distributions of charm and anticharm 
quark scattering with quarks/antiquarks (solid line) and 
gluons (dashed line) as a function of the invariant en¬ 
ergy y/s in 0-10 % central Au+Au collisions. The total 
number of charm and anticharm quark scatterings with 
light quarks and light antiquarks is about 135 and that 
with gluons is about 83. We recall that about 19 pairs 
of charm and anticharm quarks are produced in this cen¬ 
trality range; each charm or anticharm quark thus experi¬ 
ences on average 6 elastic scatterings with partons before 
it is hadronized. Since the dressed quark/antiquark mass 
is smaller than the dressed gluon mass, the peak of the 
distribution for charm and quark scattering is located at 
a lower energy compared to that of charm and gluon scat¬ 
tering. We additionally note that the number of charm 
quark scatterings with gluons is not small compared to 
that with quarks although the number of gluons is sig¬ 
nificantly smaller than that of quarks due to their larger 
mass. This is because the cross section for the charm and 
gluon elastic scattering is about twice that for charm and 
light quark elastic scattering. One characteristic of the 
partonic scattering of charm quarks is that the scatter¬ 
ing distribution in y/s has a long tail up to high invariant 
energies beyond that expected in a thermal equilibrium. 












FIG. 6: (Color online) Spatial diffusion coefficient of heavy 
quark, D a , as a function of T for p q = 0. The black solid line 
below T = 180 MeV is the hadronic diffusion coefficient [32|, 
and the red solid line above T c ~ 160 MeV our D s computa¬ 
tion in the partonic environment. The lattice calculations are 
from Ref. (49]. 



FIG. 7: (Color online) The distributions of charm and an- 
ticharm quark scattering by quarks/antiquarks (solid line) 
and gluons (dashed line) as a function of the invariant en¬ 
ergy yfs in 0-10 % central Au+Au collisions. 


This is attributed to the cross section for charm and par- 
ton elastic scattering which does not decrease with in¬ 
creasing y/s as shown in Fig. [4] Therefore, the partonic 
scattering is effective for the energy loss of charm and 
anticharm quarks at high transverse momentum. 

Fig. [8] compares the transverse-momentum spectrum 
of charm and anticharm quarks at the initial produc¬ 
tion with that at hadronization in the midrapdity interval 
(|y| < 1). We can see that the scattering with massive 


FIG. 8: (Color online) Transverse momentum spectra of 
charm and anticharm quarks at the initial production (dashed 
red line) and at hadronization (solid red line) in 0-10 % cen¬ 
tral Au+Au collisions, compared with those of thermal charm 
quarks with (dot-dashed blue line) and without (dotted blue 
line) considering the transverse flow of the bulk particles at 
T=158 MeV which is the critical temperature for the QGP 
phase transition in PHSD. 


partons softens the transverse momentum spectrum of 
charm quarks substantially. These spectra are also com¬ 
pared in Fig. [8] with the thermal spectrum of charm and 
anticharm quarks at T c with the number of charm and 
anticharm quarks being the same as in the PHSD. The 
dot-dashed blue line in Fig. [8] shows the spectrum when 
we take into account the transverse flow velocity vt of 
the bulk particles which is about 0.44 at T c in 0-10 % 
central collisions in the PHSD simulations: 


dN 

2npTdpTdy 


m T Io(oi)Ki(P), 


(7) 


where m-r = y/m% + p|,, a = pr sinh p/T and (3 = 
rriT coshp/T with p = tanh -1 vt, and K\ and Iq are 
modified Bessel functions [ho) . 

The dotted blue line shows the spectrum without in¬ 
cluding the transverse flow (v t = 0). Since the trans¬ 
verse flow velocity of charm quarks in relativistic heavy- 
ion collisions is smaller than that of the bulk particles, 
even though they are completely thermalized j5lj |. the 
transverse-momentum spectrum of thermalized charm 
quarks is between the dotted blue and the dot-dashed 
blue lines. This figure suggests that charm and anticharm 
quarks are close to thermal equilibrium at low transverse 
momentum (p? < 2 GeV) after partonic scattering while 
they are still off-equilibrium at higher transverse mo¬ 
menta since the solid red and the dot-dashed blue line 
start to deviate substantially for higher pt. 
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V. HADRONIZATION OF CHARM QUARKS 


Since the hot and dense matter created by a relativistic 
heavy-ion collision expands with time, the energy density 
of the matter decreases and the deconfined degrees-of- 
freedom hadronize to color neutral hadronic states. Once 
the local energy density in PHSD becomes lower than 0.5 
GeV/fm 3 , the partons are hadronized. First we look for 
all combinations of a charm quark and light antiquarks 
or of an anticharm quark and a light quark and calculate 
the probability for each combination to form a D or D* 
(or D s , D *) meson. The probability for a quark and an 
anti-quark to form a meson is given by 


, n 8 9m 

f(P, kp) = —- exp 


_li_ _ l+A 2 


where gM is the degeneracy of the meson M , and 

1 , , r- m 2 ki - mik 2 

p = —t={v i-r 2 ), k p = V2 




mi + m 2 


( 8 ) 


(9) 


with rrii, r * and k, being the mass, position and momen¬ 
tum of the quark or antiquark i. respectively. The width 
parameter 6 is related to the root-mean-square radius of 
the meson produced through 


(r 2 ) 


3 m\ + m\ 

2 (mi + m 2 ) 2 


( 10 ) 


and thus determined by experiment (if available). 

The degeneracy factors, gM, for D and D* mesons 
are 1 and 3, respectively. We also take into account 
higher excited states of D mesons from the Particle Data 
Group [12: Z^(2400)°, D 1 (2420)°, and D^(2460)°’ ± . 
They dominantly decay into D + tt or D* + 7r and each 
branching ratio is not known yet. Therefore, we assume 
that the higher excited state decays immediately after 
hadronization into D + ir or D* + 7r with the branch¬ 
ing ratio of 1 to 3, keeping its flavor. We further simplify 
that higher excited states have the same mass, 2460 MeV, 
since D2(2460) o,± are most abundantly produced in the 
coalescence model due to its large spin, and they are as¬ 
sumed to have the same radius as D and D* mesons. 

In the actual simulations of the PHSD we perform sev¬ 
eral tens of events in parallel to obtain smooth local en¬ 
ergy and particle number densities as a function of time; 
the charm or anticharm quark in the hadronization is al¬ 
lowed to take its partner from other events. Therefore, 
the coalescence probability given in Eq. © is divided 
by the number of parallel events such that the results 
become independent on the number of parallel runs. 

Collecting all possible combinations of a charm or an 
anticharm quark with light quarks or antiquarks and cal¬ 
culating the coalescence probability for each combination 
from Eq. ©, we obtain the probability for the charm 
or the anticharm quark to hadronize by coalescence in 
the actual space-time volume AtAxAyAz. Whether the 
charm or the anticharm quark is actually hadronized by 


coalescence is decided by Monte Carlo. Once the charm 
or the anticharm quark is decided to be hadronized by co¬ 
alescence, then we find its partner again by Monte Carlo 
on the basis on the probability of each combination in the 
selected local ensemble. In case the charm or anticharm 
quark is decided not to hadronize by coalescence, it is 
hadronized by the fragmentation method as in p+p col¬ 
lisions (see Sec. m■ Since the hadronization by coales¬ 
cence is absent in p+p collisions, it can be interpreted 
as a nuclear matter effect on the hadronization of charm 
and anticharm quarks. 




FIG. 9: (Color online) The probabilities for charm or an¬ 
ticharm quarks to hadronize to D or D meson through coa¬ 
lescence as a function of transverse momentum in 0-10, 10-40, 
40-80 % central Au+Au collisions at ^/snn =200 GeV for the 
D meson radius of 0.5 fm (upper) and of 0.9 fm (lower). 


Fig.[9]shows the probabilities for charm and anticharm 
quarks to hadronize to D and D mesons through coales¬ 
cence as a function of transverse momentum px in 0-10, 
10-40, 40-80 % central Au+Au collisions at ^/snn =200 
GeV for the D meson radius of 0.5 fm (upper) and 0.9 fm 
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(lower). We can see that the probability decreases as the 
transverse momentum of the charm or anticharm quark 
increases. The reason is that the charm or anticharm 
quark with larger transverse momentum has less chance 
to find a neighboring coalescence partner in phase space. 

Fig. [9] also shows that coalescence probability depends 
on centrality and the radius of the D meson. Since there 
are more light partons as a coalescence partner for the 
charm or anticharm quark in central collisions, the coales¬ 
cence probability increases with decreasing centrality. In 
extremely peripheral collisions, however, the coalescence 
probability becomes tiny - as in p+p collisions - and the 
nuclear modification factor of D mesons naturally ap¬ 
proaches to 1. Since the hadronization of charm quarks 
through fragmentation is not properly suited at low 
transverse momentum, several studies have forced the 
coalescence probability to be 1 at pr =0 mmmm- 
Unless only a single centrality is concerned, it is techni¬ 
cally difficult to do so for all centralities by using the same 
coalescence parameters. Moreover, if the coalescence 
probability is 1 in peripheral collisions, it will induce a 
large nuclear modification factor, although small nuclear 
matter effect is expected in such collisions. Therefore, 
we allow the fragmentation of charm quarks at low trans¬ 
verse momentum as in p+p collisions in competition with 
coalescence (for different radii of the D-mesons). In Fig. [9] 
we see that the coalescence probability is larger for the 
D meson radius of 0.9 fm than of 0.5 fm. We note that 
while the smaller radius is more physical for D-mesons, 
a somewhat larger radius has been used in several stud¬ 
ies [13, [EH to get a larger coalescence probability. 

The energy-momentum difference between the charm 
or anticharm quark and the fragmented D or D meson 
in fragmentation and the energy difference in coalescence 
are distributed equally to the surrounding partons (in the 
same cell) which are not hadronized yet to conserve the 
total energy-momentum of the system. 


VI. D MESON SCATTERING IN THE HADRON 

GAS 

The D and D* mesons produced through coalescence 
or fragmentation interact with the surrounding hadrons 
in PHSD. The presence of several resonances close to 
threshold energies with dominant decay modes involv¬ 
ing open-charm mesons and light hadrons suggests that 
the scattering cross sections of a D/D* off a meson 
or baryon, highly abundant in the post-hadronization 
medium, could manifest a non-trivial energy, isospin 
and flavor dependence. An example of these states is 
the broad scalar resonance Do (2400), which decays into 
the pseudoscalar ground state D by emitting a pion in 
the s-wave (similarly to the heavy-quark spin partner 


Ah (2420), decaying into D*^ 1 ). Moreover, the similar¬ 
ity between the A(1405) and the A c (2595) has driven the 
attention to the fact that the latter could be playing the 
role of a sub-threshold resonance in the DN system, con¬ 
nected to the latter by coupled-channel dynamics. 

All these features have been addressed within sev¬ 
eral recent approaches based on hadronic effective mod¬ 
els which incorporate chiral symmetry breaking in the 
light sector. The additional freedom stemming from the 
coupling to heavy-flavored mesons is constrained by im¬ 
posing heavy-quark spin symmetry (HQSS) |32l l+ll 1591 . 
Whereas chiral symmetry fully determines the scattering 
amplitudes of Goldstone bosons with other hadrons at 
leading order in a model independent way, by means of 
HQSS the dynamics of the pseudoscalar and the vector 
mesons containing heavy quarks can be connected, since 
all kinds of spin interactions are suppressed in the limit 
of infinite quark masses. 

Following [EiJ (and references therein), the Lagrangian 
density describing the interaction between the spin-zero 
and spin-one D mesons with the light pseudoscalar 
(Goldstone) bosons from the octet (7r, A', A, rf) reads 
C = £lo + £nlo> where the subscripts LO and NLO 
refer to the leading and next-to-leading orders in the chi¬ 
ral perturbative expansion, whereas one keeps at leading 
order in the heavy-quark expansion. The LO Lagrangian 
reads 

Clo = {WDV ll D')-m%(DD')-{WD w V lt Dy) 

+ m 2 D (D*^D/}) + ig{D *»- Du^D/}) 

+ ^(D; Ua V 0 D*J -V 0 D;u a D*J)e^ ,( 11 ) 

where D = (D°,D+,D+) and A>* = (D*°, D*+, £>*+)„ 
are the SU(3) antitriplets of spin-zero and spin-one 
D mesons, respectively, with mass mo in the chi¬ 
ral limit, and the brackets denote the trace in flavor 
space. Two interaction terms are present in the LO 
Lagrangian, containing the coupling of the D and D* 
meson fields to the light-meson axial vector current, 
Un = i (v/d^u — ud^u'), which describe, for instance, 
the decay of the D* meson into a Dtt pair. We note 
that the two terms are in principle independent but 
HQSS demands that both are related by a unique cou¬ 
pling constant g. The covariant derivative contains the 
coupling to the light-meson vector current and reads 
= 0 m “ I ( m1 A u + ud^v/) , where the Goldstone 
bosons are introduced within the non-linear realization 


1 According to quark model predictions, the excitation spectrum 

of the D / D* system consists of a triplet and a singlet of positive 

parity, with J p = 0"*", 1"*”, 2"*" and 1"*", respectively (correspond¬ 
ing to 2S+1 Lj = 3 Pj and 1 Pi). The l - * - states mix in such a way 
that the meson with the lower mass, D i(2420), becomes narrow 
and decouples from the D*7r low-energy dynamics. The higher 
one, Di(2420), on the contrary, is very broad and decays pre¬ 
dominantly in D* 7T. The tensor state, ^2(2460), is also narrow 
and its influence can be, in principle, disregarded. 
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of chiral symmetry in exponential parameterization, U = 
u 2 = exp , with 



/ 

tsJ M 

o 

+ 

s 

7T + 

I<+ 

$ = 


7r 

-A— tt 0 + 1 „ 

K° 


V 

K~ 

K° 

~7e T1 


and / the meson decay constant in the chiral limit (/ = 
93 MeV). 

The NLO Lagrangian, which we omit here for simplic¬ 
ity, introduces twelve low-energy constants (LECs), hi 
and hi (i = 0,..., 5), which have to be fixed by sym¬ 
metry arguments and known phenomenology. Imposing 
HQSS at LO reduces the number of free parameters to 
six, since hi = hi, whereas by large- N c considerations it 
can be argued that only the odd LECs contribute. 

The interaction of D mesons with baryons has been 
recently studied in different approaches which incorpo¬ 
rate heavy flavor in the meson baryon interaction, ei¬ 
ther by implementing f-channel vector-meson exchange 
mechanisms between pseudoscalar mesons and baryons, 
by extending the Jiilich meson-exchange model, or rely¬ 
ing on the hidden gauge formalism (see [so] and refer¬ 
ences therein). Among these we recourse to the model 
described in Refs. [57H591 |Hj, which generalizes the 
Weinberg-Tomozawa form of the meson-baryon interac¬ 
tion, fixed by chiral symmetry at leading order, beyond 
flavor SU(3) structure. In addition, in this model HQSS 
is fulfilled exactly whenever charm quarks participate in 
the interaction. In all charm sectors, and in particular 
in C = 1 (C = charm number), the theory accounts for 
point-like s-wave interactions between charmed mesons 
and light mesons in the pseudoscalar and vector octets, 
and the lowest lying baryons from the J p = l/2 + (nu¬ 
cleon) octet and the 3/2 + (Delta) decuplet. It reduces 
to the WT interaction prescribed by chiral symmetry 
as long as the pseudoscalar mesons (Goldstone bosons) 
are involved. The latter is accomplished by enhancing 
the standard SU(3) flavor symmetry to a SU(6) spin- 
flavor symmetry [SU(8) when charm is also considered], 
which has allowed to identify unambiguously many bary- 
onic states as dynamically generated resonances of the 
meson-baryon interaction in the light sector f57i - [59l . [61 1]. 
HQSS imposes that arbitrary rotations of the heavy- 
quark spin should leave dynamics unchanged, which pre¬ 
vents charm flavor to be exchanged between the inter¬ 
acting meson-baryon pair, naturally leading to the sup¬ 
pression of charm-exchange processes 2 . The tree-level 


2 This feature arises in other models based on meson exchange 
mechanisms due to the heavy mass of the charmed field and the 
large heavy-meson decay constants; note, however, that HQSS is 
not exactly fulfilled in most of these models fvl . 


amplitudes of the theory are given by 


V IJSC 

v ij 


(Vs) 


use 2 Vs — Mi — Mj 
U ij 


tfifi 


Ei + Mi / Ej + Mj 


2 M % 


2 Mi 


(13) 


where Mi and Ei denote the mass and CM energy of the 
baryon in channel i, respectively, /* stands for the corre¬ 
sponding meson decay constant, and jj IJSC is a matrix of 
coefficients in the coupled-channel space for given isospin, 
spin, strangeness and charm numbers. Breaking of the 
SU(8) symmetry is taken into account by using physical 
hadron masses and meson decay constants in Eq. m 
The model is not complete in the analytic sense since no 
t- or u- channel mechanisms are explored; it is, however, 
a minimal extension of the SU(3) WT interaction with 
no additional free parameters which simultaneously im¬ 
plements the relevant symmetries of QCD in the light- 
and heavy-quark sectors. 

The tree-level scattering amplitudes obtained from the 
previous approaches to the scattering of D/D* mesons 
with light mesons and baryons are unitarized along the 
right-hand cut by solving the set of coupled-channel (on- 
shell) Bethe-Salpeter equations, T = T + VGT. Here 
T (V) stands for the unitarized (tree-level) amplitude 
with matrix element T{V)lj SC , where i(j) labels the 
incident (outgoing) heavy-light meson or heavy-meson 
baryon state. G represents the two-particle propagator or 
loop function, which is regularized by subtractions or in 
dimensional regularization. The few parameters uncon¬ 
strained in the heavy-light meson Lagrangian are fixed 
by the mass difference between the D and D s mesons, 
and by the position and width of the Dq( 2400) reso¬ 
nance, which is dynamically generated by the unitarized 
interaction in the 7=1/2 Dir channel. We recall that 
the model for the meson-baryon interaction is parameter 
free, modulo regularization. The sub-threshold A c (2595) 
resonance is dynamically generated by the model, with a 
strong coupling to the DN and D*N channels in I = 0 
and a finite width from the E c 7r channel, which is open 
for decay. Several additional states in the C = 1 and 
S = 0 sector are generated by the interaction, some of 
which have been identified with experimental observa¬ 
tions whereas others are genuine predictions that could 
be detected in the forthcoming CBM and PANDA exper¬ 
iments at GSI/FAIR [58l ]. 

A selection of the D and D* cross sections off light 
mesons, nucleons and A baryons is shown in Fig. [TO] 
for several physical (charge) states. The cross sections 
exhibit a remarkable isospin dependence, as expected, 
leading to rather different shapes (energy dependence). 
For example, the Dp (2400) resonance is clearly visible 
in the D°ir + elastic reaction, which has a predominant 
7 = 1/2 component. Conversely, the D°ir~ elastic reac¬ 
tion is pure 7 = 3/2 and the cross section is not resonant 
in this case, whereas other processes exhibit deviations 
from the resonance profile due to the interference between 
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FIG. 10: (Color online) Several examples for the scattering cross sections of D and D* mesons with light mesons and baryons. 


the / = 1/2,3/2 amplitudes. We note the dip caused by 
the opening of the Dp channel, which mixes with Dir 
only at NLO in the chiral Lagrangian. The DK cross 
sections present a characteristic monotonic fall due to the 
presence of the narrow D*(2317), which lies right below 
threshold energy (compare, e.g, with the cross section of 
the pure isovector D°K ~, which is essentially flat). The 
isoscalar part of the DK amplitude contributes with a 
broad resonant state around 2700 MeV. Significantly, the 
vector partner Z?* 1 (2700) has been observed with a de¬ 
cay width of about 120 MeV. Similar cross sections are 
obtained for the D* interaction as the underlying dynam¬ 
ics is identical at LO in the heavy-quark expansion, and 
the only differences are introduced by the slightly higher 
mass of the D* (thus thresholds and resonances appear 
shifted by this amount). In addition to the subthreshold 
A c and £ c states generated from the D and D* interac¬ 
tions with baryons, the D^N and D^A cross sections 
are populated by several resonances at higher energies, 


with a remarkable isospin and coupled-channel depen¬ 
dence. In both cases the reactions with a strong isovec¬ 
tor component seem to dominate over other processes. 
We also account explicitly for the cross sections of reac¬ 
tions with antinucleons, not shown in Fig. 1101 which are 
obtained from the DN scattering amplitudes by charge 
conjugation. Finally, the cross sections for the scattering 
of D and D* with the vector mesons from the octet (e.g. 
Dp), out of the scope of the present model, have been set 
to an estimated value of 10 mb and independent of the 
collision energy. 

Fig. [Ill shows the distribution of D and D* meson scat¬ 
terings with pions and nucleons as a function of the in¬ 
variant scattering energy ^/s in 0-10 % Au+Au central 
collisions. The total number of D and D* meson scat¬ 
terings with pions is about 49 and that with nucleons or 
antinucleons is ~ 6. The numbers of D and D* meson 
scatterings with all mesons and all baryons are, respec¬ 
tively, ~ 56 and ~ 10. Considering that the number of D 
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FIG. 11: (Color online) The distributions of D and D* meson 
scattering with pions and nucleons as a function of the invari¬ 
ant scattering energy \/s in 0-10 % Au+Au central collisions. 
See the legend for the individual channels. 

and D* mesons in 0-10 % central collisions is about 30, 
excluding D s , D*, A c , and A c , each D or D* meson expe¬ 
riences on average two scatterings with a hadron until it 
freezes out. Compared to Fig. [71 the number of scatter¬ 
ings decreases rapidly with increasing scattering energy. 
This is attributed to the decrease of hadronic cross sec¬ 
tions beyond the resonance region, as shown in Fig. 1101 
Accordingly, hadronic interactions become ineffective for 
the energy loss of D and D* mesons at high transverse 
momentum. 


VII. RESULTS 




The nuclear modification of D mesons is expressed in 
term of the ratio Raa which is defined as 


R-aa(pt) = 


dN Au+Au /dp T 

N^XdN^/dpr’ 


(14) 


where V^ u+Au and N^ u+Au are, respectively, the num¬ 
ber of D mesons produced in Au+Au collisions and that 
in p+p collisions, and N A ^ AAu is the number of binary 
nucleon-nucleon collisions in Au+Au collision for the cen¬ 
trality class considered. If the matter produced in rela¬ 
tivistic heavy-ion collisions does not modify the D meson 
production and propagation, the numerator of Eq. (ED 
should be equal to the denominator. Therefore, an Raa 
smaller or larger than one implies that the nuclear matter 
suppresses or enhances D mesons, respectively. 

In order to see the effect of partonic scattering, all 
hadronic scatterings of D mesons are switched off in 
Fig. El The dotted, dashed, solid, and dot-dashed lines 
in the figure show the Raa of D° mesons in 0-10 % central 
Au+Au collisions with the cross sections for the partonic 


FIG. 12: (Color online) The Raa of D° mesons without 
hadronic scattering in 0-10 % central Au+Au collisions at 
^/snn =200 GeV for a radius of the D meson of 0.5 fm (up¬ 
per) and of 0.9 fm (lower), which are compared with the ex¬ 
perimental data from the STAR Collaboration §]. The cross 
sections for charm-parton scattering are multiplied by factors 
of 0, 0.5, 1, and 2 for the dotted, dashed, solid, and dot-dashed 
lines, respectively. 


scattering of charm being artificially multiplied by factors 
of 0, 0.5, 1, and 2, respectively. We stress that these mul¬ 
tiplication factors are introduced in order to explore the 
impact of the partonic interaction strength on the shape 
of Raa(pt )■ Note, however, that the default cross sec¬ 
tion (multiplication factor 1) is determined consistently 
on the basis of the DQPM couplings and propagators and 
thus has no free parameters. The fact that the Raa{pt ) 
is best described by the consistent cross sections (blue 
solid lines) points towards an experimental support of 
our approach. The black dotted line, where both par¬ 
tonic and hadronic scatterings are absent, shows that 
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FIG. 13: (Color online) The Raa of D° mesons in Au+Au 
collisions at ^/snn =200 GeV including partonic and hadronic 
scatterings for the D meson radius of 0.5 fm (dashed lines), 
of 0.9 fm (solid lines), and without coalescence, i.e. frag¬ 
mentation only (dotted lines), which are compared with the 
experimental data from the STAR Collaboration [§]. 


Raa approches 1 at high p T as it ought to be. However, 
the coalescence of charm quarks creates a peak in Raa 
around px = 2.5 GeV/c. This peak does not show up in 
the PHSD calculations when discarding the dissolution 
of strings to partons. In this case the Raa is unity. The 
peak emerges because the charm hadron gains transverse 
momentum in the coalescence which is absent in p+p col¬ 
lisions. The red dashed, blue solid, and green dot-dashed 
lines show that with increasing cross section for charm 
and parton scattering the charm quark looses more en¬ 
ergy at high pt■ It is also seen that the energy loss of a 
charm or anticharm quark at high pt can be dominantly 
attributed to the interaction with partons in the QGP. 
On the other hand, the Raa at low px increases with 


increasing scattering cross section. The reason is that 
a larger scattering cross section produces charm quarks 
closer to their thermal equilibrium distribution as shown 
in Fig. 1 which is enhanced at low px relative to the 
initial px distribution. 



FIG. 14: (Color online) The Raa of D° mesons including 
partonic scattering with (dashed and solid lines) and without 
hadronic scattering (dotted and dot-dashed lines) in Au+Au 
collisions at ^/snn =200 GeV for a D meson radius of 0.5 fm 
and of 0.9 fm. The experimental data are from the STAR 
Collaboration 0] . 

In Fig. Iwe show the effect of the coalescence prob¬ 
ability on the Raa of D° mesons. The dashed and solid 
lines are, respectively, the Raa for a D meson radius of 
0.5 fm and of 0.9 fm as shown in Fig. [9] For the dotted 
lines, all charm and anticharm quarks are hadronized by 
fragmentation. It is seen that for the D meson radius 
of 0.9 fm more charm quarks are hadronized through co¬ 
alescence and as a result the peak of Raa is shifted to 
higher px■ On the contrary, when all charm quarks are 
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FIG. 15: (Color online) The elliptic flow V 2 of D° mesons 
including partonic scattering with (dashed and solid lines) and 
without hadronic scattering (dotted and dot-dashed lines) in 
Au+Au collisions at ^/snn =200 GeV for a D meson radius 
of 0.5 fm and of 0.9 fm. The experimental data are from the 
STAR Collaboration 00. 

hadronized through fragmentation, the peak is shifted to 
lower p t ■ 

Finally, Fig. [TH shows the Raa and Fig. [15] the elliptic 
flow V 2 of D° mesons with and without hadronic scat¬ 
tering. We can see that the hadronic scattering plays an 
important role both in Raa and the clliptc flow V 2 ■ It 
shifts the peak of Raa to higher transverse momentum 
especially in central collisions and enhances the elliptic 
flow of final D mesons. 


VIII. SUMMARY 

We have studied charm production in relativistic 
heavy-ion collisions by using the Parton-Hadron-String 
Dynamics (PHSD) approach [27] where the initial charm 
quark pairs are produced in binary nucleon-nucleon colli¬ 
sions from the PYTHIA event generator [43| taking into 
account the smearing of the collision energy due to the 
Fermi motion of nucleons in the initial nuclei. The pro¬ 
duced charm and anticharm quarks interact with the 
dressed quarks and gluons in the QGP which are de¬ 
scribed by the Dynamical Quasi-Particle Model [38] in 
PHSD. The interactions of the charm quarks with the 
QGP partons have been evaluated with the DQPM prop¬ 
agators and couplings consistently . Furthermore, 
when extracting the spatial diffusion coefficient D s from 
our cross sections (Fig. [6]) as a function of the tempera¬ 
ture we observe a minimum of D s close to T c which is in 
line with lattice data above T c and hadronic many-body 
calculations below T c . 

We recall that the PHSD differs from conventional 
Boltzmann approaches incorporating on-shell scattering 


with pQCD cross sections in a couple of essential aspects: 

i) it incorporates dynamical quasi-particles due to the fi¬ 
nite width of the spectral functions; 

ii) it involves scalar mean-fields for the light partons that 
substantially drive the collective flow in the partonic 
phase and includes c-quark scattering with the QGP 
partons that transfers collective flow also to the charm 
quarks; 

iii) it is based on a realistic equation of state from lattice 
QCD and thus describes the speed of sound c s (T) reliably 
(without incorporating a first order phase transition); 

iv) the hadronization of ’bulk’ partons is described by 
the fusion of off-shell partons to off-shell hadronic states 
(resonances or strings) and does not violate the second 
law of thermodynamics; 

v) all conservation laws (energy-momentum, flavor cur¬ 
rents etc.) are fulfilled in the hadronization (contrary to 
some coalescence models); 

vi) the effective partonic cross sections no longer are given 
by pQCD and are evaluated within the DQPM in a con¬ 
sistent fashion and probed by transport coefficients (cor¬ 
relators) in thermodynamic equilibrium (shear- and bulk 
viscosity, electric conductivity, magnetic susceptibility, 
spatial charm diffusion coefficient etc. |35l 1401 Hi}). 

We have found that the interaction with the dy¬ 
namical partons of the QGP softens the pt spectrum 
of charm and anticharm quarks but does not lead to 
a full thermalization for transverse momenta pt > 2 
GeV/c. The charm and anticharm quarks, furthermore, 
are hadronized to D mesons either through the coales¬ 
cence with a light quark or antiquark or through the frag¬ 
mentation by emitting soft ’perturbative’ gluons. Since 
the hadronization through coalescence is absent in p+p 
collisions, it can be interpreted as a nuclear matter ef¬ 
fect on the D meson production in relativistic heavy-ion 
collisions. In the coalescence mechanism the charm or 
anticharm quark gains momentum by fusing with a light 
quark or antiquark while it looses momentum in the frag¬ 
mentation process (as in p+p reactions). This partly con¬ 
tributes to the large Raa of D mesons between 1 and 2 
GeV/c of transverse momentum. Finally, the formed D 
mesons interact with hadrons by using the cross sections 
calculated in an effective lagrangian approach with heavy 
quark spin-symmetry [32}, which is state-of-the art. We 
have found that the contribution from hadronic scatter¬ 
ing both to the Raa and to the elliptic flow of D mesons 
is appreciable, especially in central collisions, and pro¬ 
duces additional elliptic flow V 2 - 

Since the PHSD results reproduce the experimental 
data from the STAR Collaboration without radiative en¬ 
ergy loss in the pt range considered, we conclude that 
collisional energy loss is dominant at least up to pt = 
6 GeV/c in relativistic heavy-ion collisions. In our ap¬ 
proach this is essentially due to the infrared enhanced 
coupling a s (T) in the DQPM leading to large scattering 
cross sections of charm quarks with partons at temper¬ 
atures close to T c and to rather massive gluons in the 
partonic bulk matter. It will be of future interest to 
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perform a similar study at LHC energies for Pb+Pb re¬ 
actions since here a significantly larger p t range can be 
addressed and the effect of gluon bremsstrahlung might 
become important again [481 ]. Furthermore, angular cor¬ 
relations between pairs of D D mesons are expected to 
provide further valuable information. 
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